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Abstract. Applying the microcanonical definition of entropy to a weakly interacting and self- 
gravitating neutralino gas, we evaluate the change in the local entropy per particle of this gas 
between the freeze out era and present day virialized halo structures. An "entropy consistency" 
criterion emerges by comparing the obtained theoretical entropy per particle of the virialized halos 
with an empirical entropy per particle given in terms of dynamical halo variables of actual galactic 
structures. We apply this criterion to the cases when neutralinos are mostly B-inos and mostly 
Higgsinos, in conjunction with the usual "abundance" criterion requiring that present neutralino 
relic density complies with 0.2 < f2„i < 0.4 for h ~ 0.65. The joint application of both criteria 
reveals that a much better fitting occurs for the B-ino than for the Higgsino channels, so that 
the former seems to be a favored channel along the mass range of 150GeV < m ? i < 250 GeV. 
These results are consistent with neutralino annihilation patterns that emerge from recent theoretical 
analysis on cosmic ray positron excess data reported by the HEAT collaboration. The suggested 
methodology can be applied to test other annihilation channels of the neutralino, as well as other 
particle candidates of thermal WIMP gas relics . 



INTRODUCTION. 

There are strong theoretical arguments favoring lightest supersymmetric particles (LSP) 
as making up the relic gas that forms the halos of actual galactic structures. Assuming 
that R parity is conserved and that the LSP is stable, it might be an ideal candidate 
for cold dark matter (CDM), provided it is neutral and has no strong interactions. The 
most favored scenario 121 l2l Ij, yl |y, tZD considers the LSP to be the lightest neutralino 
(Xj), a mixture of supersymmetric partners of the photon, Z boson and neutral Higgs 
boson J3|]. Since neutralinos must have decoupled once they were non-relativistic, it is 
reasonable to assume that they constituted originally a Maxwell-Boltzmann (MB) gas 
in thermal equilibrium with other components of the primordial cosmic plasma. In the 
present cosmic era, such a gas is practically collision-less and is either virialized in 
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galactic and galactic cluster halos, in the process of virialization or still in the linear 
regime for superclusters and structures near the scale of homogeneity ll8Ll9L ll0ll . 

Besides the constraint due to their present abundance as main constituents of cosmic 
dark matter (Q.^\ ~ 0.3), it is still uncertain which type of annihilation cross section char- 

Xo 

acterizes these neutralinos. In this paper we present a method that discriminates between 
different cross sections, based on demanding (together with the correct abundance) that 
a theoretically estimated entropy per particle matches an empiric estimate of the same 
entropy, but constructed with dynamic variables of actual halo structures. The applica- 
tion of this "entropy consistency" criterion is straightforward because entropy is a state 
variable that can be evaluated at equilibrium states, irrespectively of how enormously 
complicated could be the evolution between each state. In this context, the two fiducial 
equilibrium states of the neutralino gas are (to a good approximation) the decoupling 
(or "freeze out") and their present state as a virialized relic gas. Considering simplified 
forms of annihilation cross sections, the joint application of the abundance and entropy- 
consistency criteria favors the neutralinos as mainly "B-inos" over neutralinos as mainly 
"higgsinos". These results are consistent with the theoretical analysis of the HEAT ex- 
periment lfllL[l2L[l3ll which aims at relating the observed positron excess in cosmic rays 
with a possible weak interaction between neutralinos and nucleons in galactic halos. 

The paper is organized as follows. In section 2 we describe the thermodynamics of 
the neutralino gas as it decouples. Section 3 applies to the post-decoupling neutralino 
gas the entropy definition of the microcanonical ensemble entropy, leading to a suitable 
theoretical estimate of the entropy per particle. In section 4 we obtain an empiric 
estimate of this entropy based on actual halo variables, while in section 5 we examine 
the consequences of demanding that these two entropies coincide. Section 6 provides a 
summary of these results. 



THE NEUTRALINO GAS 

The equation of state of a non-relativistic MB neutralino gas is 



X = 
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where m»i and are the neutralino mass and number density. Since we will deal exclu- 

Xo Xo 

sively with the lightest neutralino, we will omit henceforth the subscript «i, understand- 

Xo 

ing that all usage of the term "neutralino" and all symbols of physical and observational 
variables (i.e. Q.q, m, p, rc, etc.) will correspond to this specific particle. As long as the 
neutralino gas is in thermal equilibrium, we have 
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where g = 1 is the degeneracy factor of the neutralino species. The number density n 
satisfies the Boltzmann equation 



n + 3Hn = -(o|v|) n 2 -{n^) 



(4) 



where H is the Hubble expansion factor and (o|v|) is the annihilation cross section. Since 
the neutralino is non-relativistic as annihilation reactions "freeze out" and it decouples 
from the radiation dominated cosmic plasma, we can assume for H and (o|v|) the 
following forms 
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(o|v|) = a + b(v 2 ) : 
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where m p = 1.22 x 10 19 GeV is Planck's mass, g* = g*(T) is the sum of relativistic 
degrees of freedom, (v 2 ) is the thermal averaging of the center of mass velocity (roughly 
v 2 oc \ jx in non-relativistic conditions) and the constants a and b are determined by the 
parameters characterizing specific annihilation processes of the neutralino (s-wave or 
p-wave) [3]. The decoupling of the neutralino gas follows from the condition 



T = n(a\y\) = H, 



(7) 



leading to the freeze out temperature 7f. Reasonable approximated solutions of © 
follow by solving for Xf the implicit relation yD 
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where g^f = g*(7f) and cq ~ 1/2 yields the best fit to the numerical solution of © and 
O. From the asymptotic solution of © we obtain the present abundance of the relic 
neutralino gas [3] 
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where So ~ 4000cm 3 is the present radiation entropy density (CMB plus neutrinos), 
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1.05 x 10" 5 GeVcm- 3 . 



Since neutralino masses are expected to be in the range of tens to hundreds of GeV's 
and typically we have Xf- ~ 20 so that 7f > GeV, we can use g^f ~ 106.75 B4J] in 
equations ® - (flOb . Equation © shows how Xf has a logarithmic dependence on m, 



while theoretical considerations Jill II II EE! related to the minimal supersymetric 
extensions of the Standard Model (MSSM) yield specific forms for a and b that also 
depend on m. Inserting into (l^b- dTTjb the specific forms of a and b for each annihilation 
channel leads to a specific range of m that satisfies the "abundance" criterion based on 
current observational constraints that require 0.1 < CLq < 0.3 and h « 0.65 Un- 
suitable forms for (o\v\) can be obtained for all types of annihilation reactions 10]. If 
the neutralino is mainly pure B-ino, it will mostly annihilate into lepton pairs through 
t-channel exchange of right-handed sleptons. In this case the cross section is p-wave 
dominated and can be approximated by © with lEl[l4. 151 
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where m/ is the mass of the right-handed slepton (m/ ~ m [4]) and a 2 = gi/4% ~ 
0.01 is the fine structure coupling constant for the t/(l)y gauge interaction. If the 
neutralino is Higgsino-like, annihilating into W-boson pairs, then the cross section is 
s-wave dominated and can be approximated by © with 0,(3 [3] 
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where m w = 80.44 GeV is the mass of the W-boson and = g^/^n ~ 0.03 is the fine 
structure coupling constant for the SU (2)l gauge interaction. 

In the freeze out era the entropy per particle (in units of the Boltzmann constant k B ) 
for the neutralino gas is given by [8.ll0fl9ll 



P+P 
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where we have assumed that chemical potential is negligible and have used the equation 
of state (HJ). From ® and (fT3l . it is evident that the dependence of Sf- on m will be 
determined by the specific details of the annihilation processes through the forms of a 
and b. In particular, we will use (fTTT) and (fT2l to compute Sf from (l8l-(fT3l. 



THE MICROCANONICAL ENTROPY 

After the freeze out era, particle numbers are conserved and the neutralinos constitute 
a weakly interacting and practically collision-less self-gravitating gas. This gas is only 
gravitationally coupled to other components of the cosmic fluid. As it expands, it expe- 
riences free streaming and eventually undergoes gravitational clustering forming stable 
bound virialized structures lll0L l9L ll6Lll7ll . The evolution between a spectrum of density 
perturbations at the freeze out and the final virialized structures is extremely complex, 
involving a variety of dissipative effects characterized by collisional and collision-less 
relaxation processes iflH ITvL [Till . However, the freeze out and present day virialized 
structures roughly correspond to "initial" and "final" equilibrium states of this gas. 



Therefore, instead of dealing with the enormous complexity of the details of the in- 
termediary processes, we will deal only with quantities defined in these states with the 
help of simplifying but general physical assumptions. 

The microcanonical ensemble in the "mean field" approximation yields an entropy 
definition that is well defined for a self-gravitating gas in an intermediate scale, be- 
tween the short range and long range regimes of the gravitational potential. This inter- 
mediate scale can be associated with a region that is "sufficiently large as to contain a 
large number of particles but small enough for the gravitational potential to be treated 
as a constant" 11611 . Considering the neutralino gas in present day virialized halo struc- 
tures as a diluted, non-relativistic (nearly) ideal gas of weakly interacting particles, its 
microcanonical entropy per particle under these conditions can be given in terms of the 
volume of phase space lfl7ll 
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where V and E are local average values of volume and energy associated with the 
intermediate scale. For non-relativistic velocities v/c <C 1, we have V « \/n <=< m/p 
and E mv 2 /2 m/x. In fact, under these assumptions the definition (fl4l . evaluated at 
the freeze out, is consistent with © and (fT3t . and so it is also valid immediately after 
the freeze out era (once particle numbers are conserved). Since (fl4l is valid at both the 
initial and final states, respectively corresponding to the decoupling (sj , Xf, rif) and the 
values (s (h \ x' h \ n' h> ) associated with a suitable halo structure, the change in entropy per 
particle that follows from (fl4l between these two states is given by 
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where ([T3l can be used to eliminate Sf- in terms of x^. Considering present day halo 
structures as roughly spherical, inhomogeneous and self-gravitating gaseous systems, 
the intermediate scale of the microcanonical description is an excellent approximation 
for gas particles in a typical region of ~ 1 pc 3 within the halo core, near the symmetry 
center of the halo where the gas density enhancement is maximum but spacial gradients 
of all macroscopic quantities are negligible lfl9H2(ill . Therefore, we will consider current 
halo macroscopic variables as evaluated at the center of the halo: jcj?, ri^. 

In order to obtain a convenient theoretical estimate of from (fT5l) . we need to relate 
with present day cosmological parameters like Q.q and h. Bearing in mind that density 
perturbations at the freeze out era were very small (dnf/rif < 10 -4 , iE 111 E3]X the 
density is practically homogeneous and so we can estimate it from the conservation of 

particle numbers: rif = no (1 +Zf) 3 , and of photon entropy: g^fSf = g*oSo (1 +Zf) 3 , valid 
from the freeze out era to the present for the unperturbed homogeneous background. 
Eliminating (1 +Zf) 3 from these conservation laws yields 
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where g*o = g*(r o CMB ) ~ 3.91 and r o CMB = 2.7 K. Since for present day conditions no/n 
po/p? and p = 

Pcrit ^oh 2 , we collect the results from (O and write (JT3I) as 
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Therefore, given m and a specific form of (o|v|) associated with a and Z?, equation (fTTl) 
provides a theoretical estimate of the entropy per particle of the neutralino halo gas that 
depends on the initial state given by Xf- in © and (fT3l . on observable cosmological 
parameters h and on generic state variables associated to the halo structure. 



THEORETICAL AND EMPIRIC ENTROPIES 



If the neutralino gas in present halo structures strictly satisfies MB statistics, the entropy 
per particle, s®, in terms of p£' = rnnj 1 and x£> = mc 2 /{k B T^), follows from the well 



known Sackur-Tetrode entropy formula (21J 
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Such a MB gas in equilibrium is equivalent to an isothermal halo if we identify ||2: 



— = — — = a 2 , 



(19) 



where a 2 is the velocity dispersion (a constant for isothermal halos). 

However, an exactly isothermal halo is not a realistic model, since its total mass 
diverges and it allows for infinite particle velocities (theoretically accessible in the 
velocity range of the MB distribution). More realistic halo models follow from "energy 
truncated" (ET) distribution functions I23I I24I l28ll that assume a maximal "cut 

off" velocity (an escape velocity). Therefore, we can provide a convenient empirical 
estimate of the halo entropy, s'^\ from the microcanonical entropy definition (fl4l in terms 
of phase space volume, but restricting this volume to the actual range of velocities (i.e. 
momenta) accessible to the central particles, that is up to a maximal escape velocity 
v e (0). From theoretical studies of dynamical and thermodynamical stability associated 
with ET distribution functions 1231 l24l I25L l28l I27L l29l 12611 a nd from observational data 
for elliptic and LSB galaxies and clusters ll30Ll3llll9ll32l l33h. it is reasonable to assume 



v 2 (0) = 2 1^(0) I ~ aa 2 (0), 12<a<18, 
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where <J>(r) is the newtonian gravitational potential. We have then 
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where we used jcj? = c 2 /o^(0) as in ( fT9l) . As expected, the scalings of d2Tt are identical 
to those of (fT8l . Similar entropy expressions for elliptic galaxies have been examined in 
& 

Comparison between obtained from (OTt and from (fTTt leads to the constraint 
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which does not depend on the halo variables xf, p*', hence it can be interpreted as 
the constraint on = 5/2+Xf- that follows from the condition sf\ th = s^l^. Since 
we can use © and ([TUl) to eliminate h 2 Q.Q, the constraint (l2"2l) becomes a relation 
involving only x^, m, a, b, a. This constraint is independent of ©, which is another 
(independent) expression for Sf = 5/2+Xf, but an expression that follows only from 
the neutralino annihilation processes. Therefore, the comparison between s'^\ ih and s^\ em , 
leading to a comparison of two independent expressions for Sf, is not trivial but leads 
to an "entropy consistency" criterion that can be tested on suitable desired values of 
m, a, b, a. This implies that a given dark matter particle candidate, characterized by m 
and by specific annihilation channels given by Xf- through ©, will pass or fail to pass this 
consistency test independently of the details one assumes regarding the present day dark 
halo structure. This is so, whether we conduct the consistency test by comparing © and 
(l2"2l or (fTTt and (OTT) . However, the actual values of for a given halo structure, whether 
obtained from (I2TT) or from (fTTt . do depend on the precise values of p*' and xf. Since the 
matching of either (HJ) and (l22l or (fTTt and d2TT) shows a weak logarithmic dependence 
on m, the fulfillment of the "entropy consistency" criterion identifies a specific mass 
range for each dark matter particle. This allows us to discriminate, in favor or against, 
suggested dark matter particle candidates and/or annihilation channels by verifying if 
the standard abundance criterion © is simultaneously satisfied for this range of masses. 



TESTING THE ENTROPY CONSISTENT CRITERION 

Since we can write (|2"2"1) as: 



ln(/i 2 Q ) =7.57-Jc f + ln ((Uf) 
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this constraint becomes a new estimate of the cosmological parameters h?Q.Q, given 
as in terms of a structural parameter of galactic dark matter halos, a, the mass of the 




FIGURE 1. Figures (a) and (b) respectively correspond to the Higgsino and B-ino channels. The 
shaded regions display Q.q vs m from our entropy criterion d23l with the solid curve giving Q.q from the 
cosmological abundance criterion (|9jl, in all cases for h — 0.65. The horizontal dashed lines give current 
estimates of Q.q = 0.3 ±0.1. It is evident that only the B-ino channels allow for a simultaneous fitting of 
both the abundance and the entropy criteria. 



neutralino, m, and the temperature of the neutralino gas at freeze out, Xf-. This last 
quantity depends explicitly not only on m, but also on its interaction cross section, and 
hence on the details of its phenomenological physics viz ©■ 

At this point we consider values for the constants a and b that define the interaction 
cross section of the neutralino, and use (1231 to plot Q.q as a function of m in GeV's. 
Using h = 0.65 and given the uncertainty range of a, we will obtain not a curve, 
but a region in the Qrj — m plane. Considering first condition (fT2"l) . corresponding to 
Higgsino-like neutralinos, leads to the shaded region in figure la. On this figure we 
have also plotted the relation which the abundance criterion © yields on this same 
plane. Firstly, we notice that the mass range that results from our entropy criterion 
intersects the one resulting from the abundance criterion. However, it is evident that 
within the observationally determined range of Q.q (the horizontal dashed lines 0.2- 
0.4), there is no intersection between the shaded region and the abundance criterion 
curve. This implies that both criteria are mutually inconsistent, thus the possibility that 
Higgsino-like neutralinos make up both the cosmological dark matter and galactic dark 
matter appears unlikely. 

Repeating the same procedure for mainly B-ino neutralinos, (fTTTl yields figure lb. In 
this case, we can see that the abundance criterion curve falls well within the shaded 
region defined by the entropy criterion. Although we can not improve on the mass 
estimate provided by the abundance criterion alone, the consistency of both criteria 
reveals the B-ino neutralino as a viable option for both the cosmological and the galactic 
dark matter. 

It is also interesting to evaluate (l2TTl and (fTTt for the two cases of neutralino channels: 
the B-ino and Higgsino, but now considering numerical estimates for jc' h) and p m that 
correspond to central regions of actual halo structures. Considering terminal velocties 
in rotation curves we have vj- erm ~ 2cj? (0) , so that x£_ ~ 2(c/vt e rm) 2 > while recent 
data from LSB galaxies and clusters H32L 133. 20, 36] suggest the range of values 
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FIGURE 2. Figures (a) and (b) respectively correspond to the Higgsino and B-ino channels. The 
figures display 4 h, | cm from d20i — d2 li (gray strip), s*'| [h from dl7> for h = 0.65 and the uncertainty strip 
ilo = 0.3 ±0.1 (thick curves) and 4 h, | MB from Jl 8I > (crosses), all of them as functions of log 10 m. The 
vertical strip marks the range of values of m that follow from d^l-jlOt for the same values of i2o and h. 
It is evident that only the B-ino channels allow for a simultaneous fitting of both the abundance and the 
entropy criteria. 



0.01M Q /pc 3 < p£> < lM /pc 3 . Hence, we will use in the comparison of (fT7t and (l2lT) 
the following numerical values: p' c hl = 0.01M Q /pc 3 = 0.416 GeV/cm 3 and x V = 2 x 10 6 , 



typical values for a large elliptical or spiral galaxy with vterm — 300 km/sec ||35|,|20, 3 



Figure 2a displays s^\ a and s'p\ cm as functions of log 10 m, for the halo structure described 
above, for the case of a neutralino that is mostly Higgsino. The shaded region marks 
sp\ m given by d2TT) for the range of values of a, while the vertical lines correspond to the 
range of masses selected by the abundance criterion © for = 0.2, 0.3, 0.4. The solid 
curves are s^\ a given by ( fT71) for the same values of Q.q, intersecting the shaded region 
associated with (OTT) at some range of masses. However, the ranges of coincidence of 
a fixed ( fT7t curve with the shaded region (OTT) occurs at masses which correspond to 
values of Q.q that are different from those used in <TTvT> . that is, the vertical lines and solid 
curves with same Q.q intersect out of the shaded region. Hence, this annihilation channel 
does not seem to be favored. 

Figure 2b depicts the same variables as figure 2a, for the same halo structure, but 
for the case of a neutralino that is mostly B-ino. In this case, the joint application of 
the abundance and entropy criteria yield a consistent mass range of 150 GeV < m«i < 
250 GeV), which allows us to favor this annihilation channel as a plausible dark matter 
candidate, with m lying in the narrow ranges given by this figure for any chosen value of 



CIq. As noted above, the results of figures la and lb are totally insensitive to the values 
of halo variables, xf® and p™, used in evaluating (ETT) and (flTt . Different values of these 
variables (say, for a different halo structure) would only result in a relabeling of the 
values of sf along the vertical axis of the figures. 



CONCLUSIONS 

We have presented a robust consistency criterion that can be verified for any annihilation 
channel of a given dark matter candidate proposed as the constituent particle of the 
present galactic dark matter halos. Since we require that the empirical estimate sj? | m 
of present dark matter haloes must match the theoretical value s^\, derived from the 
microcanonical definition and from freeze out conditions for the candidate particle, the 
criterion is of a very general applicability, as it is largely insensitive to the details of 
the structure formation scenario assumed. Further, the details of the present day halo 
structure enter only through an integral feature of the dark halos, the central escape 
velocity, thus our results are also insensitive to the fine details concerning the central 
density and the various models describing the structure of dark matter halos. A crucial 
feature of this criterion is its direct dependence on the physical details (i.e. annihilation 
channels and mass) of any particle candidate. 

Recent theoretical work by E. A. Baltz et al. JTlll confirmed that neutralino anni- 
hilation in the galactic halo can produce enough positrons to make up for the excess 
of cosmic ray positrons experimentally detected by the HEAT collaboration lfl2l [Till . 
Baltz et al. concluded that for a boost factor B s ~ 30 the neutralinos must be primarily 
B-inos with mass around 160 GeV. For a boost factor 30 < B s < 100, the gaugino- 
dominated SUSY models complying with all constraints yield neutralino masses in the 
range of 150 GeV < m^i < 400 GeV. On the other hand, Higgsino dominated neutrali- 
nos are possible but only for B s ~ 1000 with masses larger than 2 TeV The results that 
we have presented in this paper are in agreement with the predictions that follow from 



I llll . as we obtain roughly the same mass range for the B-ino dominated case (see figure 
lb) and the Higgsino channel is shown to be less favored in the mass range lower than 
TeV's. 

We have examined the specific case of the lightest neutralino for the mostly B-ino 
and mostly Higgsino channels. The joint application of the "entropy consistency" and 
the usual abundance criteria clearly shows that the B-ino channel is favored over the 
Higgsino. This result can be helpful in enhancing the study of the parameter space of 
annihilation channels of LSP's in MSSM models, as the latter only use equations ® 
and d9b — <TT0b in order to find out which parametersyield relic gas abundances that are 
compatible with observational constraints SSSHIaED- However, equations © and 
(l^l)- (fT(l by themselves are insufficient to discriminate between annihilation channels. 
A more efficient study of the parameter space of MSSM can be achieved by the joint 
usage of the two criteria, for example, by considering more general cross section terms 
(see for example [3]) than the simplified approximated forms (fTTIl and (fT2l . This work 
is currently in progress. 
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